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Several kinds of spin density wave (SDW) states with both quarter-filled band and dimerization 
are reexamined for a one-dimensional system with on-site, nearest-neighbor and next-nearest- 
neighbor repulsive interactions, which has been investigated by Kobayashi et al. (J. Phys. Soc. 
Jpn 67 (1998) 1098). Within the mean-field theory, the ground state and the response to the 
density variation are calculated in terms of phase variables, and <f>, where 9 expresses the charge 
fluctuation of SDW and <j> describes the relative motion between density wave with up spin and 
that with down spin respectively. It is shown that the exotic state of coexistence of 2/cf-SDW 
and 2&F-charge density wave (CDW) is followed by 4/cf-SDW but not by 4/cf-CDW where &F 
denotes a Fermi wave vector. The harmonic potential with respect to the variation of and/or 
<f> disappears for the interactions, which lead to the boundary between the pure 2/cf-SDW state 
and the corresponding coexistent state. 

KEYWORDS: spin density wave, quarter-filled band, charge density, phase variable, long range Coulomb interaction 
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§1. Introduction 

Spin density wave (SDW) state has been studied 
extensively for quasi-one-dimensional organic conduc- 
tors, (TMTSF) 2 X and (TMTTF) 2 X, which are known 
as Bechgaard saltHfEP Since the energy band of the 
salt is anisotropic and quarter-filled, the SDW state is 
determined essentially by the nesting along the one- 
dimensional direction where the corresponding Fermi 
wave vector, fc F , is fc F — 7r/4 in the unit of lattice con- 
stant. The ground state of such a SDW, which has been 
studied by NMR experiment ,t!lt$ exhibits the periodic ar- 
ray of magnetic moment along the one-dimensional axis 
where the periodicity is four times as large as the lat- 
tice constant and the magnitude varies every two lattice 
sites. The degree of dimerization, which reduces the sys- 
tem effective Ly to the half- filled one, has been estimated 
in these salt.EP Another noticeable feature of this SDW 
state is the charge fluctuation, which is caused by the 
deformation of the phase of the density wave, as seen in 
the nonlinear electric conductivity.!]) 

Recently the X-ray scattering investigation has shown 
that 2/cp-SDW coexists with 2fc F -charge density wave 
(CDW) in the SDW state of (TMTSF) 2 PF 6 at temper- 
atures just below the onset temperature of the SDW 
ordering.EP The 4fc F satellite reflection is also observed 
for both (TMTSF) 2 PF 6 and (TMTTF) 2 Br. Further the 
coexistent state of 2fc F -SDW, 2fc F -CDW and 4fc F -CDW 
has been verified in the temperature range 4-12 K corre- 
sponding to SDW1 phase while CDW is absent in SDW3 
phases at lower temperatures!? 

Theoretical study for the coexistence of SDW and 
CDW has been explored_based on the mean-field the- 
ory. Seo and FukuyamaEP have found the coexistence 
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of 2fc F -SDW and 4fc F -CDW by examining the extended 
Hubbard model, which takes into account both nearest- 
neighbor repulsive interaction and dimerization. The 
SDW state, which is followed by 4fc F -CDW, appears 
when the magnitude of nearest-neighbor interaction is 
larger than a critical value. Based on this model, 
the collective mode describing the charge fluctuation 
was calculated. The excitation spectrum, which has 
a gap due to commensurability, becomes gapless for 
the nearest-neighbor interaction leading to the onset of 
the coexistence. EP The noticeable state of the coexis- 
tence of 2fc F -SDW and 2fc F -CDW has been demonstrated 
by Kobayashi et al. _svho introduced the next-nearest- 
neighbor interaction.EP They examined the SDW state 
as the function of both the nearest-neighbor interaction 
(V a ) and the next-nearest-neighbor interaction (V" 2 ) and 
found the phase diagram, in which such a coexistent state 
is located in the region with large V a and large V 2 . Thus 
the model with the long range Coulomb interaction re- 
veals the ground state with variety of SDW states. How- 
ever it is not yet clear if 4/c F -SDW state and/or 4fc F - 
CDW state can coexist in addition to the 2/c F -CDW. 

In the present paper, such a calculation-is performed 
for a model proposed by Kobayashi et al.f3 within the 
mean-field theory. The ground state is examined in terms 
of two kinds of phase variables 9 and 4>, where 6 is a phase 
of 2fc F -SDW and (f> is the difference between a phase 
of density wave with up spin and that with down spin. 
Further, as a first step for studying the density fluctua- 
tion, we examine the increase of energy, which is caused 
by varying these two kinds of phases from those of the 
ground state. In §2, formulation is given for the model, 
which comprises dimerization and alternating transfer 
energy (t a and tb) with repulsive interaction (U), the al- 
ternating nearest-neighbor interaction (V a and Vb) with 
the fixed Vb/V a and the next-nearest-neighbor interac- 
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tion (V2). The order parameters is examined in terms 
of phase variable 9 and <j>. In §3, the ground state of 
the SDW is calculated to obtain a phase diagram on the 
plane of V 2 /V a and V a /tt,. Based on the study of the 
ground state, the property of the density variation is in- 
vestigated by calculating the energy as a function of 9 
and <f>. Discussion is given in §4. 

§2. Formulation 

We study the SDW state with the quarter-filled band 
by use of a one-dimensional extended Hubbard model 
with dimerization and long range Coulomb interaction. 
The model is represented as 

N 

# = - E E (* - 1- 1 )^) { C l C 3+i,° + h-c.) 

+ H iat , (2.1) 

N N 

= uJ2 n jmi + E ( y ~ (-1)^) n 3 n i+l 



3 = 1 



3 = 1 



N 



+ V lJ2 U 3 n 3+2 
3 = 1 



(2.2) 



where t a = t + t d , t b = t — t d , V a = V + SV and 
V b = V - SV. The quantity C) a {C ja ) denotes the 
creation (annihilation) operator of the electron at the 
j-th lattice site with spin a. The quantity N is the 
total number of the lattice site. n 3<T — Gj a C ]a and 
n j = n iT + n jl- We take t, fee and the lattice con- 
stant as unity. The alternation of the transfer integrals 
given by t ± td comes from dimerization, which induces 
the energy, t^- Coupling constants for the Coulomb in- 
teraction are defined by U, V±SV and V 2 , which denote 
the on-site repulsive interaction, the nearest-neighbor in- 
teraction and the next-nearest-neighbor interaction, re- 
spectively The quantity SV also originates in the dimer- 
ization. By noting that the electron band is at quarter- 
filling, order parameters of density wave can be written 
as (m = 0, 1, 2, 3), 
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(2.3a) 

C L C k+rnQ ^ MF , (2.3b) 



where Qq = 2fep = ir/2 with fcp being the Fermi wave 
vector. The spatial dependence of electron densities with 
up and down spin, which has the periodicity being four 
times as large as the lattice constant, is obtained from 
the Fourier transform of eqs. (2.3a) and (2.3b). Then 
we examine eqs. (2.3a) and (2.3b) under the condition 
that So = 0, D = 1/2, S Qa = S* 3Qq , D Qo = D* Qq , 
S 2Q„ = S 2Qa = S 2 and D 2Qo = D* Qa = D 2 . In terms 
of the mean-field (MF) given by eqs. (2.3a) and (2.3b), 
Hamiltonian (eq. (2.1)) is expressed as 
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where e k = —2t cos k. Since the nesting vector is given by 
Qo = 7r /2, the eigenvalue of eq. (2.4) exhibits four energy 
bands E n (k), (n = 1,2,3,4), where E^k) < E 2 (k) < 
E 3 (k) < En(k). In the ground state, the lowest energy 
band, E\{k), is fully occupied and others are empty due 
to the quarter-filling. By solving eqs. (2.3a) and (2.3b) 
self-consistently, the ground state energy per site, E g , is 
obtained as 

^ = ^E E ^( fc ) 

a 0<k<Q o 
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where the Brillouin zone is reduced to < k < Q - In 
terms of eq. (2.5), we define the quantity, SE, as 



5E = E„ — -Enormal 



(2.6) 



which denotes the difference between the ground state 
energy of the ordered state and that of the normal state 
i.e., Sj = Dj ■ = ( j = 1,2). In eq. (2.6), £ nOT mai ( 
= -2y/2t/n + U/16 + V/A + V 2 /4) is the energy of the 
normal state. 

Now we examine the SDW state from the aspect of 
the phase variable. Since Sq = S^q o and Dq — D^q q , 
eqs. (2.3a) and (2.3b) are expressed as 



D ie w 



(2.7a) 
(2.7b) 



where Si( > 0) and Di( > 0) are real. Quantities 9 
and 6' are phase factors. Within the present numerical 
calculation, we find always Sq /Si = iDq /Di for the 
ground state. Thus we calculate the energy under such 
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a restriction, i.e., eq. (2.7b) is replaced by 

D Qo = Ae i(e ~ V2) ■ (2.8) 

In terms of eqs. (2.3a), (2.3b), (2.7a) and (2.8), the 
electron density for spin a at the j'-th site is expressed 
as 

(n ja ) = l/A + sgn(a)S 1 cos(Q r j +9)+D 1 sm(Q r j +9) 



+ -(D 2 + sgn(a)S 2 ) cos(2Q rj) 



(2.9) 



The respective electron density is rewritten as 
(n n ) = 1/4 + Asm(Q r J + 9 + <f>) + Bcos(2Q r,) , 

(2.10a) 

{n n ) = 1/4 + Asm(Q rj + 9 -(f)) + Ccos{2Q Q r j ) , 

(2.10b) 



where A = y/S'f + D'{, <f> = tan-^Si/Di) , B = 
(D 2 + S 2 )/2 and C = {D 2 -S 2 )/2. Equations (2.10a) and 
(2.10b) imply that the density wave with up (down) spin 
has a phase 9 + <p (9 — <fi). The quantity, 0, expresses 
the relative motion between the density wave with up 
spin and that with the down spin. The idea, that <f> 7^ 
leads to the coexistence of 2fc F -SDW and 2fc F -CDW, has 
been already asserted by OverhauserO The quantity, 
9, denotes a phase representing the charge density fluc- 
tuation.. pyhiriL. gives rise to the sliding motion of the 

SDW.BBBBO) 

The calculation of the energy can be implemented by 
two kinds of methods, i.e., parameters determined self- 
consistently are given either by Si, S 2 , D±, D 2 and 9 or 
by A, S 2 , D 2 , 9 and 4>. For the simplicity, the former 
choice is taken for the calculation of the ground state 
while the latter one is used for studying the effect of the 
variation of phases 9 and <p around the ground state. 
We define V as the energy in the presence of density 
waves with S n Q and D n Q (n=l,2 and 3), which are 
introduced as the trial function. From eqs. (2.4), (2.7a) 
and (2.8) together with (j> = tan _1 (S , 1 /Z? 1 ) , V is written 
as 

V- 1 (H * U V 

V — — \-"MF 



N x ' 16 4 
), is performed by eq. 



(2.11) 
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(2.4), in which 



The average, 
quantities S nl ^ 
Note that eqs. (2.3a) and (2.3b) are equivalent to the 
condition that 



and D n n are retained as parameters. 



dV _ dV_ 

dS 3 ~ dDj ~ ' 

^ = 
89 ' 



(.7 = 1,2) , (2.12) 



(2.13) 



and TV 1 V . (rij a ) = 1/4. Then we introduce a quantity, 
V(9, <fi), which is defined by 



V{9,(j>) = V\ 



dV _ OV _ dV _ n 



(2.14) 



The minimum of V{9, <fi) is given by 9 = 9q and <fi = <f>Q, 
i.e., quantities f?o and 4>q are phases for the ground state. 
Thus, we can evaluate the increase of energy by V(9, <fi) 



when the phase variables are varied from the ground 
state. Especially the limiting behavior for small vari- 
ation of these phases is evaluated by expanding V(9, <j>) 
in terms of 9 — 9q(= rcos tp) and <p ~ <Po{= rsin-0). For 
r -C 1, one obtains 

V(0,4>) =V(6 ,<j> ) + C(V)r a + ». , (2.15) 

where C(ip) denotes the coefficient of harmonic potential. 

§3. SDW States and Phase Variables 

3.1 Ground state 

We examine the possible ground state, which gives the 
lowest energy in the solution of the self-consistency equa- 
tions. From the numerical calculation, we found three 
kinds of solutions for SDW states shown in Fig. 1: pure 
2-fc F -SDW (I), a coexistent state of 2fc F -SDW and 4fc F - 
CDW (II) and that of 2fc F -SDW, 2fc F -CDW and 4fc F - 
SDW (III). The pure 2Jc F -SDW state is obtained for the 
interaction of only ll£3 A coexistent state of 2A: F -SDW 
and 4A; F -CDWEP appears for V lager than a critical value. 
An exotic state of the coexistence of 2& F -SDW apd 2fc F - 
CDW has been obtained for large V and large V 2 S3 This 
state differs from the state III in the sense that 4fc F -SDW 
also coexists in the state III. It seems to be reasonable 
that the state III could be possible if the order parameter 
for 4fc F -SDW is treated in a perturbational method . 

In the present calculation, we examine self-consistency 
equations by varying V a /tb and V 2 /V a . Parameters aie 
chosen as t a /t b = 1.1, V b /V a = 0.8 and U/t b = 4,Ef 
where the magnitude of dimerizjation seems to be rea- 
sonable for the Bechgaard salts B Figure 2(a) displays 
the V^/ifc-dependence of Si, S 2 and D\, which denote 
amplitudes for 2fc F -SDW, 4/fc F -SDW and 2fc F -CDW re- 



(I) 



1 


I 1 


I 










1 

* 


' T 



(II) 



—t 1 











(III) 



Fig. 1. The electron density < Tij a > as the function of the lattice 
site Tj for three kinds of SDW states: pure 2-fcp-SDW (I), a 
coexistence of 2fcp-SDW and 4fcp-CDW (II) and a coexistence of 
2fc F -SDW, 2fc F -CDW and 4fc F -SDW (III). The variation satisfies 
< rtj+4 >=< rij > where the solid (dashed) arrow denotes the 
magnitude of electron density with up (down) spin. 
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Fig. 2. (a) The V a /t^-dependence of order parameters Si, S2 and 
D\ for V^/Va = 0.75. In the inset, order parameters, Si (solid 
curve) and D2 (dashed curve) arc shown for V^/Va = and 
0.3. Order parameters are zero if not shown explicitly, (b) The 
V a /^-dependence of 6 (solid curve) for Vi/tb = and 0.3 and 
that of <f> (dashed curve) for Vi/tb = 0.75. 



spectively . For V 2 /V a = 0.75, the state I moves to the 
state III at V a /tb = 1.66, where both D\ and S 2 starts 
increasing continuously with a slight enhancement of Si . 
Such a simultaneous appearance of D\ and S 2 originates 
in the fact that the term proportional to D1S2S1 exists in 
the expansion of the Ginzburg-Landau free energy. Note 
that D 2 = 0, i.e., 4^-CDW is absent for V 2 /V a = 0.75. 
In the inset, the V^/tb-dependence of order parameters 
for = and 0.3 is shown where solid curve and 

dashed curve denote S\ and D 2 respectively. The tran- 
sition from the state I into the state II takes place at 
V a /t b = 0.76 and 1.13 for V 2 /V a = and 0.3 respec- 
tively where S 2 — Di = for any V a jtb- It is found 
that Si in the state I does not depend on both V a and 



Fig. 3. The V^/Va-dependence of order parameters, Si, S2, D\ 
and D2, for V a /tt =1.6. 



V 2 and that the state II is followed by 4fcF-CDW. In Fig. 
2(b), quantity 9 corresponding to V 2 /h = and 0.3 is 
shown by the solid curve. When the state I moves to the 
state II, 9 decreases rapidly from 7r/4 and <f> remains at 
4> = tt/2 for both states. The dashed curve denotes (j> 
for V 2 /tb = 0.75, which shows a transition from the state 
I into the state III. The quantity <f> decreases from 7r/2 
while 9 — tt/A for both states I and III. Note that the 
coexistence of 2/cf-SDW and 2fcp-CDW in the state III 
originates in </> with < <f> < tt/2. 

Figure 3 shows the successive transition from the state 
II, the state I and the state III with increasing V 2 /V a - 
Order parameters Si, S 2 , Di and D 2 are shown as a 
function of V 2 /V a with the fixed V 2 /tb =1.6. A transition 
from the state II into state I occurs at V 2 /V a = 0.48 while 
a transition from the state I into state III occurs at V 2 /V a 
= 0.78. The ^-dependence of D 2 , D\ and S2 indicates 
the transition being of the second order. The quantities 
Di and S2 appear at the same V 2 /V a , as seen in Fig. 2. 
A linear dependence of the variation of Si , which is found 
just below the transition at V 2 /V a = 0.48, comes from 
the effect of higher order with respect to D 2 . In the state 
I, Si does not depend on V 2 since the pure 2&F-SDW is 
determined only by the [/-term, i.e., the density-density 
interaction of both V ± SV-term and V2-term does not 
give rise to the coupling between the order parameter of 
pure 2fcp-SDW and others. The appearance of Di and 
D 2 is understood as follows. The distance for the next- 
nearest-neighbor interaction, V 2 , corresponds to the 2fcp- 
periodicity. Then the 2/cp-CDW is induced for V 2 larger 
than a critical value, e.g., Di appears for V 2 /t>1.3 and 
irrespective of the magnitude of V. In a similar way, one 
can understand the role of the nearest-neighbor interac- 
tion, V, which induces D 2 for large V a (or V) due to 
the 4fcp-periodicity. The interval region in Fig. 3 corre- 
sponding to the state I decreases with increasing V and 
disappears for V (and/or V a (= V + SV) ) larger than a 
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Fig. 4. Phase diagram on the plane of V2/V a and V a /t b with 
t a /h = 1.1, Vt/Va = 0.8 and U/t b = 4. The closed circle given 
by (V2,V a ) = (V2 lC ,V a ,c) denotes an intersection of the solid 
curve and the dotted curve. The solid curve (dotted curve), 
which shows the boundary between the region I and the region 

II (the region I and the region III) becomes the actual boundary 
for V a < V a ,c- The boundary between the state II and the state 

III is given by the dashed curve. The thick dash-dotted curve 
corresponds to the dotted curve with S2 = ut3 and the thin 
dash-dotted curve is explained in the inset of Fig. 7(b). 



critical value. 

In Fig. 4, a phase diagram of SDW states is shown 
on the plane of Vi/V a and V a /tb for U/tb = 4, which is 
equal to U/t= 3.81. The closed circle corresponding to 
the intersection of the solid curve and the dotted curve 
denotes a critical point given by (V2, c /V a , V a . c /tb)(— 
(0.576,2.077), at which all the energies of states I, II 
and III become equal. The solid curve (the dotted curve) 
shows the boundary between the sate I and the state II 
(the state I and the state III). The dashed curve is the 
boundary between the state II and the state III. On the 
dashed curve, the energy of the state II becomes equal 
to that of the state III, while amplitudes of both order 
parameters are finite. For V a < V a , c , the second order 
transition occurs on both the solid curve and the dotted 
curve while the first order transition from the state II 
into the state III is obtained on the dashed curve. We 
note that the state III pith the assumption of S2 = has 
been already obtainedtlP where the resultant boundary 
between the state I and the state III is shown by the 
thick dash-dotted curve in Fig. 4. Compared with the 
present boundary of the dotted curve, the region for the 
state III is suppressed in the absence of 6*2. The thin 
dash-dotted curve, which is not a boundary, is explained 
in the discussion of the inset of Fig. 7(b). 

The first order transition shown by the dashed curve 
of Fig. 4 is examined by calculating the energy gain in 
the ordered state, 5E, which is defined by eq. (2.6). In 
Fig. 5, the V2/V a -dependence of SE is shown for V a /tb — 
1.6 and 3. The solid curve, the dotted curve and the 



Fig. 5. The V^/Va-dependence of energy <5E ( eq. (2.6)) with the 
fixed Va/tb = 1-6 and 3, for the state I (solid curve), the state 
II (dotted curve) and the state III (dash-dotted curve). 



dash-dotted curve denote 5E for the state I, II and III 
respectively. The quantity SE for the state I, which is 
independent of the magnitude of V a /tb, is largest com- 
pared with those of the state II and the state III. The 
case of V a jtb = 1.6 corresponds to Fig. 3, which shows a 
successive transition given by the state II, the state I and 
the state III with increasing V^/V^. For V a /tb = 3, one 
obtains a first order transition from the state II into the 
state III at V2/V a — 0.45, which is shown by the arrow. 
Based on these calculations, we find the V2/T4 and V a /tb 
dependences of 8E in Fig. 4 as follows. The quantity 
5E in the region I does not depend on both V2/V1 and 
V a /tb while SE in region II (III) decreases motonically 
with moving away from the boundary of the solid curve 
(the dotted curve). 

3.2 Property around the ground state 

In the previous section, the ground sate has been stud- 
ied by calculating order parameters where their phases 
9(— * #0 ) an d 4>{-^ 4>o) are also determined so as to mini- 
mize the energy. We examine, in this section, the prop- 
erty around the ground state by calculating V(9, 0) given 
by eq. (2.14) which denotes the energy obtained by the 
phase variation. The study is focused on the state around 
the boundary given by the solid curve and the dotted 
curve in Fig. 4. 

Figure 6(a) shows the ^-dependence of V(6, 7r/2) for 
three cases which are located close to the boundary of 
the solid curve of Fig. 4. Solid curve, dotted curve 
and dash-dotted curve are calculated for V a /tb = 0.5 
(1), 1.13 (2) and 1.5 (3) respectively with the fixed 
Vi/V a = 0.3 where curve (2) corresponds to the bound- 
ary between the state I and the state II. The quantity 
d 2 V(9 : n/2)/d9 2 at 9 = 9 , which denotes the second 
derivative of V(9,ir/2) with respect to 9 at the ground 
state, is positive for curves (1) and (3) but becomes zero 
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&H/(&P/4) V 2 /V a 
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Fig. 6. (a) The ^-dependence of V(0,ir/ 2 ) for ^V*6 = °- 5 Wi 
1.13 (2) and 1.5 (3) with the fixed V 2 /V a = 0.3. (b) The in- 
dependence of V(7r/4,<£) for Vo/tj, = 1.4 (1), 1.66 (2) and 1.8 
(3) with the fixed V 2 /V a = 0.75. 



Fig. 7. (a)The V 2 /Va-dependence of C(0) (eq. (2.15)) for 
V /t& = 1, 2 and 3. (b)The Vi/Va-dependence of C(vr/2) for 
Va/tb = 1, 2 and 3 where the right vertical axis is for Va/% = 1. 
In the inset, the details of Vg/Va -dependence of C'(n/2) is shown 
for Va/tb = 1 and 1.6. 



for curve (2). Such a-result has been already found for 
V 2 = and SV = oH From Fig. 6(a), V(6,w/2) can 
be written as V(9,tt/2) ~ const. + W c cos 40 + Wd sin 20, 
where 1F C depends on both V a and V2 and Wa is mainly 
determined by f d . With increasing V a (V2), W c de- 
creases and becomes negative (increases and approaches 
to the limiting value with V a = 0). In Fig. 6(b), the 
(^-dependence of V (ir / 4, <fi) is shown for the three cases 
near the boundary of dotted curve of Fig. 4. Solid curve, 
dotted curve and dash-dotted curve are calculated for pa- 
rameters with (V a /t b ,V 2 /V a ) = (1.4,0.75) (1), (1.66,0.75) 
(2) and (1.8,0.75) (3), which correspond to the state I, 
the boundary and the state III respectively. For curve 
(2), the second derivative of V(7r/4, cj>) with respect to <j> 
becomes zero at <f> — ir/2 while harmonic potential exists 
for curve (1) and (3). Thus it turns out that the har- 



monic potential with respect to <f> vanishes only at the 
boundary between the state I and the state III. 

It is of interest to examine how the harmonic potential 
for 9 ((f)) reduces to zero when parameters vary around 
the boundary given by the solid curve (the dotted curve) 
in Fig. 4. We calculate C(ip) given by eq. (2.15), which 
denotes a coefficient of harmonic potential with respect 
to 9 and </) around 9o and (f>o, i.e., those of the ground 
state. Quantities C(0) and C(tt/2) correspond to coeffi- 
cients for 9 and <j> respectively. In Fig. 7(a), the ^/In- 
dependence of C(0) is shown for V a /tt = 1, 2 and 3 
respectively. For V a /tb = 1, C(0) becomes zero at VijV a 
= 0.22, which corresponds to the boundary between the 
state I and the sate II. Such a behavior is similar to the 
y-dependence of the charge gap, which has been shown 
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for V 2 = 0.0) For V a /t b = 2, (7(0) decreases to zero at 
V 2 /V a = 0.56, i.e., the boundary between the state I and 
the state II and exhibits a cusp at V 2 /V a = 0.61, i.e., 
the boundary between the state I and the state III. A 
shoulder is found in C(0) for larger V 2 /V a . For V a /t b — 
3, C(0) exhibits a jump, at which a first order transition 
takes place due to the boundary between the state II and 
the state III. When V 2 increases further, C(0) takes a 
minimum and increases rapidly. In Fig. 7(b), the V2 /In- 
dependence of C(tt/2) is shown for the same parame- 
ters as Fig. 7(a). For V a /t b — 1, C(n/2) takes a cusp 
at V 2 /V a = 0.22 corresponding to a transition from the 
state II to the state I. With increasing V 2 /V a , C{tt/2) de- 
creases rapidly and becomes zero at V 2 /V a = 0.90 where 
the quadratic behavior is seen around zero. The fact 
that such a zero value does not imply the boundary is 
explained at the end of this paragraph. For V a /tb = 2. 
C{tt/2) shows a small cusp at V 2 /V a = 0.56 correspond- 
ing to a transition from the state II to the state I. A 
zero of C(tt/2) at V 2 /V a — 0.61 comes from the bound- 
ary between the state I and the state III. For V a /t b = 
3, C(n/2) takes a jump at V 2 /V a — 0.45, which origi- 
nates in the first order transition between the state II 
and the state III. In the inset, we show the detail be- 
havior around the zero of C(%/2) for V a /t b = 1 and 1.6. 
Quadratic behavior around the zero is seen for V a /tb=X. 
There are two kinds of zeros for Va/tb — 1.6 where the 
larger V 2 corresponds to the transition from the state I 
to the state III and exhibits a linear dependence around 
the zero. For V a /t b = 1.6, the quadratic behavior is 
also found around the smaller V 2l at which the actual 
transition does not occur. Such a fact is understood as 
follows. The effect of small variation of <p on the ex- 
pectation value of eq.(2.11) is essentially determined by 

D Qo = N- 1 E<r=T4 Y.-K<k<ir ( C L C k+Q ,v) , where 
the average <> is performed over eq.(2.4) with S m Q 
and D m Q replaced by trial fields, e.g., D\ oc cos</) and 
Si oc sin 4>. The ground state of the state I is the pure 
2fcp-SDW with <f) — (f)Q — (ir/2) and then quantities Dq 
and Dq q in the second line of eq. (2.4) vanishes due 
to D\ oc cos0 and eq. (2.7b). When <fi is varied from 
4>o = (7r/2), the magnitude of Dq and Dq q becomes 
finite and the resultant Dq becomes finite generally. 
However there is a case where Dq = even for Dq 7^ 
due to vanishing of the coefficient, e.g., U/2 — 2V2 = 
for SV = 0. Such a case is examined explicitly in the 
next section. The location for C(ir/2) = correspond- 
ing to Dq = is shown by the thin dash-dotted curve 
in Fig. 4, where the vanishing of C(ir/2) exists only in 
the region I. The trajectory (thin dash-dotted curve) is 
nearly parallel to the dotted curve and crosses with the 
solid curve at (V a /t b ,V 2 /V a ) ^ (1.6,0.5). 

Further we examine C(ip) with arbitrary ip, which in- 
duces the variation of both 9 and (p. In Fig. 8, C(tp) 
is calculated for {V a /t b ,V 2 /V a )= (0,0), (2,0), (2,0.6), 
(2,0.8) and (2,1), where C{ip + tt/2) = C(ip). The am- 
plitude of C(ip) is large for (0, 0) (shown by a vertical 
axis of right hand side ) but becomes small in the pres- 
ence of both V a and V 2 as seen from (2, 0) and (2,1). 
The minimum of C(ip) is found at ip = for (0,0), (2,0) 
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Fig. 8. The V-dependence of C(i>) for (V a /t b ,V 2 /V a )= (0,0), 
(2,0), (2,0.6), (2,0.8) and (2,1). The right vertical axis is for 
(0,0). 



and (2,1) while it is found at ip = tt/2 for (2,0.6). Note 
that the minimum for (2,0.8) is located between rp = 
and ip = tt/2, where such a continuous variation of the 
minimum is seen for 0.80 < V 2 /V a < 0.86 and V a /t b = 2 
in the region III. The minimum of C(ip) at ip = tt/2 is 
obtained for parameters (V 2 /V a ,V a /t b ) close to the dot- 
ted curve in Fig. 4. The quantity C(ip) exhibits a small 
■(/'-dependence for (2,0.8). We comment on the behav- 
ior of C(ip) in Fig. 4. Within the numerical accuracy 
of the present calculation, we find C(ip) = at a criti- 
cal point given by (V a /t b ,V 2 /V a ) = (2.077,0.576). The 
small magnitude of C(ip) shown for (2,0.6) is obtained for 
parameters of (V a /t b , V 2 /V a ), which are located near the 
critical point and also in the region enclosed by the dot- 
ted curve, solid curve and thin dash-dotted curve. Thus 
it is expected that the fluctuations for 9 and cf> become 
large around the solid curve and the dotted curve respec- 
tively and that the fluctuations for both 9 and <f> become 
large around the critical point of the intersection of these 
two curves . 

Now we examine both 9 and ^-dependences of 
V(9, (f>) around the ground state by choosing param- 
eters (V a /t b ,V 2 /V a )= (0,0), (2,0), (2,0.6), (2,0.8) and 
(2,1), which are the same as Fig. 8. In Fig. 9, the 
^-dependence of V(9, 4>o) is shown with the fixed 4> = 4>q 
where <P /tt/2 = 1 for (0,0), (2,0) and, (2,0.6) and 
0o/tt/2 = 0.537 and 0.499 for (2,0.8) and (2,1) respec- 
tively. There is a periodicity given by V(9 + tt, (pa) — 
V(0, <po). The ^-dependence of V{9, <p ) is small for (0,0) 
and (2,0.6) ( i.e., in the region I of Fig. 4 ) while it is 
large for (2,0), (2,0.8) and (2,1) (i.e., in the region of II 
and III of Fig. 4). The latter result comes from the fact 
that the variation of 9 needs a large amount of energy 
due to the presence of the charge density of D 2 (D±) in 
the region II (the region III). Figure 10 shows the corre- 
sponding (^-dependence for V(9 , <p) where 9q/{tt/4) = 1 
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Fig. 9. The 0-dependence of V(d, cpo) for (V a /*6, V 2 /V a )= (0,0), 
(2,0), (2,0.6), (2,0.8) and (2,1) where 0o/tt/ 2 = 1 for (0,0), (2,0) 
and, (2,0.6) and 4> /n/2 = 0.537 and 0.499 for (2,0.8) and (2,1) 
respectively. 
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Fig. 10. The (/--dependence of V(0 O ,4>) for (V a /t b , V 2 /V a )= (0,0), 
(2,0), (2,0.6), (2,0.8) and (2,1) where 8 /n/4 = 1 for (0,0), 
(2,0.6), (2,0.8) and (2,1) and e/yr/4 = 0.035 for (2,0). 



for (0,0), (2,0.6), (2,0.8) and (2,1) and 0/(tt/4) = 0.035 
for (2,0) where V(e ,<t>) = V(6 , tt/2 - 4>) = V(6 a , + tt). 
In the state III, 0o differs from tt/2, which causes a small 
increase of V(9o, 4>) for the variation of </>. The variation 
of V(6q,<P) is large for (2,0) and (0,0) and is small for 
(2,0.6), (2,0.8) and (2,1). This indicates a fact that the 
fluctuation of <fr in the state I and II is small and that in 
the state III is large. The magnitudes of V(7r/4, <j>) for 
(0,0), (2,0.6), (2,0.8) and (2,1) becomes equal at (j) = tt/2, 
since all these states reduce to the state I with the energy 
being independent of both V and V-i. 
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Fig. 11. (a)Phasc diagram on the plane of V 2 /V and V/t where 
the solid curve ( the dotted curve ) denotes the boundary be- 
tween the state I and the state II ( the state I and the state 
III). Parameters are taken as (t d /t,SV/V) = (0,0) (1), (0.1,0) 
(2), (0,0.1) (3), (0.1,0.1) (4). (b)Phase diagram on the plane of 
V-2.lt and V/t where U/t = 4 and (t d /t,SV/V) = (0.1,0.1). The 
notation of curves are the same as Fig. 4. In the inset, another 
phase diagram is shown for U/t = 4 and (tj, 5V) = (0,0). 



§4. Discussion 

In the present paper, we have examined several kinds 
of SDW states in a one-dimensional quarter-filled elec- 
tron system with long range Coulomb interaction within 
the mean-field theory. The phase diagram of the SDW 
state with the states I, II and III were obtained on the 
plane of Vi/V a and V a /tb- The variation of the energy 
from that of the ground state was calculated in order to 
examine the role of phase variables of the SDW states. 

Here we examine the effect of SV and t d on the bound- 
aries in the phase diagram of Fig. 4. In Fig. 11(a), the 
several results of boundaries are shown where the solid 
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curve ( the dotted curve ) denote the boundary between 
the state I and the state II ( the state I and the state III) . 
Parameters are taken as (i d /i, SV/V) = (0,0) (1), (0.1,0) 
(2), (0,0.1) (3), (0.1,0.1) (4) where curve (4) corresponds 
to Fig. 4. Since solid curve (2) ( solid curve (1)) is the 
the same as solid curve (4) ( solid curve (3)), the bound- 
ary between the state I and the state II depends on id but 
does not depend on SV. The energy gain by the forma- 
tion of 4fcp-CDW decreases by the dimerization but not 
by the alternation of nearest-neighbor interaction. From 
the dotted curves where curves (3) and (4) are located 
below curves (1) and (2), it is found that the energy gain 
by the formation of 2fcp-CDW is increased by SV. The 
dotted curve moves slightly upward in the presence of 
id- We note that the comparison of solid curve (4) and 
solid curve (3) (dotted curve (4) and dotted curve (2)) 
demonstrates explicitly the argument by Kobayashi et 
al&3 that the region II is enlarged by in the absence of 
id ( and that SV increases the region III). 

In order to understand the coexistent state clearly, an- 
other phase diagram is calculated on the plane of V 2 /t 
and V/t where V = (V a + V b )/2 and t = (t„ + t b )/2. 
by choosing another parameters for both axes. In Fig. 
11(b), a phase diagram corresponding to Fig. 4 is shown 
with U/t = 4 and (t d /t,SV/V) = (0.1,0.1). The point, 
which shows the intersection of the solid curve and the 
dotted curve, is given by (V 2 /t,V/t)=(l.l5,2.0). We find, 
within the numerical accuracy, that the boundary be- 
tween the state I and II (the solid curve) is expressed 
as V = V% + const.. Such a ^-dependence is under- 
stood from the coefficient of D 2 q in eq. (2.4), i.e., 
(U/2 - 2V + 2V 2 ). The quantity, V, which is needed to 
create 4fc F -CDW, is replaced by V - V 2 for V 2 ^ 0. The 
dotted curve denotes a critical value of V 2 / t, above which 
the 2/cp-CDW appears. The tangent of the dotted curve 
depends on SV. In the inset, a phase diagram is shown 
for U/t = 4 and (t d ,SV) = (0,0). The dotted curve 
becomes straight due to SV — 0, since the boundary be- 
tween the state I and the state III are determined only 
by V 2 /t in the absence of id and SV. It turns out that 
the 2fcp-CDW is induced by the next-nearest-neighbor 
interaction but not by the nearest-neighbor interaction. 
The second term of eq. (2.4) shows that the coefficients, 
(U/2 — 2V 2 ), of Dq changes the sign at V 2 /t — 1 for 
U/t = 4 and that the 2fcp-CDW is induced for larger 
V 2 /t e.g., Di ^ for V 2 /t}^1.3. Further we comment 
on the condition for C(ir/2) = 0, which occurs on the 
thin dash-dotted curve in Figs. 4 and Fig. 11(b). The 
inset of Fig. 11(b) displaying C(n/2) = at V 2 = U/A, 
is interpreted from the fact that the coefficient of Dq 
in the second line of eq. (2.4) becomes zero at V% = U /4 
resulting in the zero value of 2fcp-CDW. 

Finally we discuss the experimental results based 
on our result. In the present calculation, the state, 
whichjexhibits the coexistence of the 2/cp-SDW and 2fcp- 
CDWEJ followed by 4fc F -SDW, has been obtained for 
a model with the on site, nearest-neighbor and next- 
nearest-neighbor interactions, i.e., the purely electronic 
interactions. The experimental results of the coexistence 
of 2&F-SDW and 2fc F -CDW in TMTSF-salts-f seems to 



be consistent with the calculation. However the exper- 
iment indicates also the 4/cf-CDW, which is absent in 
the present calculation. Such a discrepancy between the 
experiment and the theory may suggest a role of the 
electron-phonon interaction, which gives rise to the den- 
sity wave in addition to the purely electronic interaction. 
Further, it will be an interesting problem to elucidate if 
4fcp-SDW obtained in the present calculation is observed 
in the coexistent state of the Bechgaard salts. 
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